EPJ manuscript No. 

(will be inserted by the editor) 



A circle swimmer at low Reynolds number 

Rodrigo Ledesma-Aguilar 1 , Hartmut Lowen 2 , and Julia M. Yeomans 1 

1 The Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 1 Keble Road, Oxford OX1 3NP, United Kingdom 

2 Institut fur Theoretische Physik II: Weiche Materie, Heinrich-Heine-Universitat Diisseldorf, Universitatsstrafie 1, D-40225 
Diisseldorf, Germany 

Received: date / Revised version: date 

Abstract. Swimming in circles occurs in a variety of situations at low Reynolds number. Here we propose 
a simple model for a swimmer that undergoes circular motion, generalising the model of a linear swimmer 
proposed by Najafi and Golestanian (Phys. Rev. E 69, 062901 (2004)). Our model consists of three solid 
spheres arranged in a triangular configuration, joined by two links of time-dependent length. For small 
strokes, we discuss the motion of the swimmer as a function of the separation angle between its links. 
We find that swimmers describe either clockwise or anticlockwise circular motion depending on the tilting 
angle in a non-trivial manner. The symmetry of the swimmer leads to a quadrupolar decay of the far flow 
held. We discuss the potential extensions and experimental realisation of our model. 

PACS. 05.40.Jc Brownian motion - 47.15. G- Low-Reynolds-number (creeping) hows - 47.61-k Micro- 
and nano- scale flow phenomena - 82.70.Dd Colloids 

1 Introduction 

The physics of microswimmers is a rapidly advancing held, 

• for recent reviews see [TJ[31[31|3] . One of the simplest models 

, of a microswimmer was proposed by Najafi and Golesta- 
nian [5] in 2004: it consists of three aligned spheres that 
are linked by rigid rods whose lengths change in time be- 
tween two values. Moving with a periodic motion which 
breaks the time-reversal symmetry, this simple swimmer 
experiences a net propulsion along the rod orientation. A 
number of other swimmer models have been proposed sub- 
sequently [6l[71l51[gi HUlllllll21ll31ll4j . most of which lead to 
propulsion along a linear trajectory. 

However, there are many examples in nature showing 

. circle swimming rather than swimming along a straight 

- line. On a planar substrate, certain bacteria [TollHTllTTlfTSl 

. [TW2"0"] and spermatozoa [2"Tll2"2"Il2"5] swim in circles. More- 
over, spherical camphors have been shown to exhibit cir- 
cular swimming when confined to an interface |24) . Re- 
cently, catalytically or thermally driven colloidal particles 
with an asymmetric shape have been prepared [25] also 
resulting in circular motion on a substrate [2"S] . 

Even though circle-swimming is frequent in the pres- 
ence of solid surfaces, curved trajectories are also very 
common in the bulk |27) . This can be attributed to asym- 
metries in the swimming stroke that can result in both 
translational and rotational modes of motion [35] . There- 
fore, while linear swimming occurs for highly symmet- 
ric swimmers, more generally swimming can occur along 
curved or circular trajectories. 



The modelling of circle swimmers is much less ad- 
vanced than that of their linearly moving counterparts. 
For instance, Dunstan et al. |29j considered a swimmer 
model of two spheres with different radii, which is a linear 
swimmer in a bulk huid but yields circle swimming close 
to a solid surface. Shum et al. [30] have implemented a 
detailed model of a flagellate that exhibits circle swim- 
ming at surfaces. For circle swimmers in the bulk, mini- 
mal rotor models .7,8,31, 32 or very coarse-grained driven 
Brownian particle models |33U34l[3"5] have been proposed. 
For the latter, the particles proceed with both an effective 
translational and angular propagation velocity and experi- 
ence additional Brownian fluctuations. The deterministic 
(noise-free) trajectory in two dimensions is a closed cir- 
cle. However, a more detailed model which resolves the 
hydrodynamic details of the swimming strokes is missing. 

In this paper we close this gap and generalise the lin- 
ear model of Najafi and Golestanian [5 to a simple circle 
swimmer. In order to do so, we consider three spheres 
joined by two links which are tilted relative to each other 
and perform the stroke as in the Najafi-Golestanian model. 
By using both analytical and numerical methods, we show 
that the resulting motion is a closed circle which depends 
on the swimmer angle, /3, which characterises the separa- 
tion between the links. We focus on the resulting radius of 
the trajectories, i?t, as a function of /3 for strokes of small 
amplitude. Interestingly, we find that swimming occurs 
predominantly along one direction of the circular trajec- 
tory (anticlockwise for the specific configuration studied 
here), except for a small range of angles, where the swim- 
mer reverses the sense of rotation. Such a behaviour is ro- 
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Fig. 1. Model of a circle swimmer. Three spherical beads of 
identical radius R are joined by two links of variable length, l\ 
and ?2- The angle between the links, /3 7^ 0, allows the swimmer 
to move in a curved trajectory. 



bust to changes in the model parameters. We further anal- 
yse the velocity field produced by this simple circle swim- 
mer and find a marked inverse-power decay at large dis- 
tances. We recover the expected quadrupolar far-field be- 
haviour [8] , with the magnitude of the velocity field decay- 
ing as the inverse cube of the distance from the swimmer. 
This is a consequence of the symmetry of the swimming 
stroke, which is invariant under combined time-reversal 
and parity transformations [371138] . For asymmetric swim- 
ming strokes, which do not posses the time-reversal and 
parity symmetry, we recover a decay consistent with a 
dipolar velocity field, as expected. 

Our simple model can serve as a starting point for fur- 
ther analytical and numerical studies. These can include 
single circle swimmers in confinement [331134] and shear 
flow 36 as well as the scattering [3 8 11313] and synchro- 
nization 40,41,42 of two circle swimmers and the (still 
unknown) collective properties of many circle swimmers, 
for example swarming and vortex formation. 

The rest of the paper is organised as follows: in Sec. [2j 
we describe and define the model. A discussion of the ge- 
ometry of the resulting trajectory is performed in Sec. [3] 
while the velocity fields in the surrounding fluid are dis- 
cussed in Sec. 2] Finally we conclude in Sec. [5] 



2 Model 

Our model swimmer consists of three beads connected by 
two massless links. These are separated by the swimmer 
angle, /3, and have lengths, l\ and Z 2 , which are known 
functions of time, and hence determine the swimming stroke 

The configuration of the swimmer can be characterised 
by the position vectors of the beads, r^, which are re- 
lated by the conditions ri 2 = l%ti, and r 2 3 = £ 2 t 2 , where 
Vij = v,j — Ti and is the unit tangent vector to link Zj. 
Since, by symmetry, motion can only occur in the plane 
defined by ti and t 2 , the orientation of the swimmer is 
determined by a single angle. Here we choose the polar 
angle associated with ti, 6q; the polar angle associated 
with t 2 hence obeys 62 = 0\ + /?■ The swimming stroke 
imposes kinematic conditions for the bead velocity vec- 
tors, Vi = Ti, where the dot indicates differentiation with 
respect to time. For the tangential motion, corresponding 



to the contraction and extension of the links, 

(v 2 - vi) • ti = ii, 



and 



h 



(v 3 - v 2 ) • t 2 
while for the angular velocities 

(v 2 -vi)-ni (v 3 -v 2 )-n 2 



(1) 
(2) 

(3) 



where h — 



by 



The bead dynamics in the overdamped limit is given 



i=i 



(4) 



where the velocity of bead i results from the force acting 
on each bead, Fj, mediated by the hydrodynamic interac- 
tion tensor, Hy, summed over all beads. For a Newtonian 
fluid of viscosity 77, the hydrodynamic interactions can be 
described using the Oseen tensor, 



H, 



1/877777^- (1 



jTij/r, 



if i = 3, 
if i ^ .1, 



(5) 



which is valid in the limit R/rij <C 1 . 

To complete the model, we impose force- free and torque- 
free conditions on the swimmer, 



and 



3 
i=l 



x F, 



0. 



(6) 



(7) 



Eqs. (HJ)-© constitute a linear system for the bead ve- 
locities and forces, whose solution depends on time through 
h(t), h(t) and /3(f) - which are prescribed functions of 
time - and through 6i(t), which reflects the dependence 
on the particular frame of reference and satisfies 9\ = 
(v 2 — Vi) • h\/l\. Once 9\(t) is known, the evolution of 
the swimmer configuration and force distribution can be 
obtained by time integrating the Vj and Fj. 

Following Refs. [5] and [8] we consider a four-step swim- 
ming stroke, where the links contract and expand alter- 
nately at constant velocity W from an initial extension D 
with a change in length e. This has the advantage of sim- 
plifying the analytics considerably, although other choices 
of the swimming stroke, consisting for instance in a contin- 
uous sinusoidal link variation, are also possible. The time 
evolution of the length of the links is chosen to be 



h{t) 



Wt 

e 

e + Wt 



if 
if 
if 
if 



< t < P/4, 
P/A<t<P/2, 
P/2<t< 3P/4, [0) 
3P/4 <t<P, 
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Fig. 2. Piece- wise swimming stroke. The links of the swimmer, 
of maximum extension D, contract and expand out of phase 
by a small length e (exaggerated in the figure) over a period 
P. 

and 

{D if < t < P/4, 

D-Wt if P/4<i<P/2, 

D-e if P/2 < i < 3P/4, W 

D-e + Wt if 3P/4 < t < P, 

and is repeated subsequently over a period P. As shown in 
Fig. [3J the swimming stroke dictated by Eqs. © and © 
breaks the time-reversal symmetry, and thus results in a 
net propulsion of the swimmer. 

3 Swimmer trajectories 

We treat the problem analytically by first solving the lin- 
ear system ([T])-([7]), we subsequently expand the Vi in pow- 
ers of R/D (up to first order), and e/D and (3 (up to 
third order) and restrict ourselves to the case where /3 
is constant in time. We are interested in the trajectory 
described by the swimmer, which can be characterised by 
angular and translational displacements over one period of 
the swimming stroke, AO and Ar, respectively. Given that 
we expect only a small angular displacement within each 
step of the stroke cycle, it is sensible to shift to the frame 
of reference where 0\ (0) = at each step and perform an 
expansion in powers of Q\ . Keeping only the leading order 
term in the expansion (constant angular velocity approx- 
imation) and adding the contribution of each step of the 
swimming stroke, we obtain, for the angular displacement, 




For the translational displacement vector of the center of 
mass, Ar = J Q P | Y^h=i v i^t = Axe x + Aye y , we obtain 
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Fig. 3. Circular trajectories for swimmers at constant swim- 
mer angle, /3. Parameter values are e/D = R/D — 10 _1 . 



and 




Eqs. (|10p - (fT2)l reduce to the result reported by Earl et 
al. [5] for vanishing /?, where the swimmer moves along a 
linear trajectory with a displacement per swimming stroke 

*"*(s)((5)" + (5)>- <-) 

For small, but finite, /3 the trajectory is no longer lin- 
ear; according to Eqs. (Tl7J|) - (rT^|) the swimmer undergoes a 
small positive rotation, and a translation along the x and 
y directions, thus describing a curved trajectory. While 
AO and Ay are odd functions of (3, Ax is an even function 
of the swimmer angle; this is consistent with the symme- 
try of the trajectory under the transformation j3 — > —fi, 
which corresponds to a reflection about the x-axis and 
has the effect reversing the sense of rotation and motion 
along the y direction while keeping the translation along 
x unchanged. 

Given that for a given f3 value AO and \Ar\ are con- 
stants, the stroke-averaged trajectories form equilateral 
chains, becoming regular polygons whenever AO = 2i:/m 
for integer m. Furthermore, due to the smallness of the 
angular displacement, the trajectories approach closed cir- 
cles with curvature 

1 AO 15 / /£\ /_£\ 2 (±\ 3 \ §_ 
~R~t ~ \M ~ 56 ^ + \d)~\d) + \d) JD' 

In order to verify this assertion, and to explore the 
full range of /3, we carry out numerical simulations for the 
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Fig. 4. (a) Normalised angular displacement as a function of the swimmer angle /3. Numerical results for the full Oseen-tensor 
hydrodynamics are carried out for e — 10~ 2 (solid line) and e = 10 _1 (dots). The dashed curve corresponds to the perturbative 
result. Inset: angular displacement during each step of the swimming stroke as a function of /3. (b) Translational displacement 
as a function of /3 for e = 10 _1 (solid line). The dashed curve corresponds to the perturbative result. 



motion of the swimmer, integrating the Oseen-level hydro- 
dynamics over time. Numerical simulations are performed 
following the algorithm proposed in Ref. [8]. At a given 
timestep, the swimmer shape is first updated according 
to the prescribed swimming stroke, while the position of 
its centre of mass and orientation are kept constant. We 
then use an iterative algorithm to find the position and 
orientation that satisfy the force-free and torque-free con- 
ditions. Swimmer trajectories, with the swimmer super- 
imposed at an arbitrary time, are shown in Fig. [3] As 
expected, trajectories are very close to circular. Fig. H] 
shows plots of AO and |Z\r| as a function of /3, where 
we also plot the perturbative result, given by Eqs. (fTU|)- 
(p~2j) . The angular displacement, depicted in Fig. EJa), 
shows a good agreement with the analytics up to angles 
as large as (3 ks 40°. As suggested by Eq. (fTU|) . results 
follow the same master curve for different deformations 
(e/D = 10 _1 and e/D = 1CP 2 ) when rescaling by the am- 
plitude (R/D)(2(e/D) 2 + 3(e/D) 3 ). This indicates that 
the shape of the curve is a function of f3 only and, con- 
sequently, that the location of maxima and minima is in- 
dependent of both R/D and e/D. Fig. HJb), shows the 
magnitude of the translational displacement, which cap- 
tures the main qualitative features of the numerical result. 
The slight discrepancy at small /3 can be attributed to cor- 
rections in R/D and e/D as discussed in Ref. [S] for linear 
swimmers. 

While their trajectories are always circular, the sense 
of rotation of the swimmers changes depending on /3, as 
shown in Fig. 0Ja). Such a dependence results from the 
competition between the four steps in the swimming stroke. 
During step 1, l\ contracts and bead 3 moves to the left, 
experiencing a drag pointing to the right (see Fig. [T]). As 
a consequence, there is a torque acting on l\ that causes a 
negative rotation of the swimmer. A similar reasoning can 



be used to conclude that the angular displacement must 
be positive for step 2. Step 3 can be mapped onto step 
2 by performing combined time-reversal and parity trans- 
formations, and therefore gives rise to the same positive 
angular displacement. Similarly, step 4 can be mapped to 
step 1. This is verified in the inset of Fig. HJa), where 
we superimpose the angular displacement at each step of 
the swimming stroke as a function of f3. Given that for 
steps 2 and 3 the relative distances between the beads are 
smaller than for steps 1 and 4, the angular displacement 
tends to be larger for the former. The net displacement 
is therefore positive for a wide range in f3. For the range 
107° < j3 < 120° the numerics show that steps 1 and 4 
dominate, causing negative net angular displacements. 

Based on their sense of rotation, we can divide swim- 
mers into three main groups. The first group, correspond- 
ing to positive AO, is delimited by two vanishing points, 
located at f3 — 0° and j3 « 107°, with a maximum located 
at /3 ~ 53°. Within the same range of angles, the trans- 
lational displacement, \Ar\, decreases with increasing j3 
to a minimum located at (3 ks 103° and then increases 
again. For the second group of swimmers the angular dis- 
placement is negative, as shown in Fig.^a), and also has a 
non-monotonic behaviour, here AO decreases for j3 > 107° 
up to a minimum value, roughly located at (3 = 113°, and 
then increases again, vanishing at j3 « 120°. In this same 
range the translational displacement increases monoton- 
ically for angles larger than 103°. For j3 > 120°, corre- 
sponding to the third group of swimmers, both the an- 
gular and translational displacements increase monoton- 
ically. These results are summarised in Fig. [5J where we 
depict the reciprocal radius of the swimmer trajectory, 
D/R t , as a function of the swimmer angle, (3, (in units 
of the rest link length, D). The sign of D/R t reflects the 
sense of motion along the trajectory, being clockwise for 
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Fig. 5. Normalised curvature of the swimmer trajectory as a 
function of the swimmer angle, j3. The solid line corresponds 
to the numerical solution of the Oseen-tensor hydrodynam- 
ics while the dashed curve corresponds to the perturbative re- 
sult. The arrows indicate the /3 values for which swimmers 
describe a linear motion. The vertical asymptote indicates a 
purely rotating swimmer. Model parameters are R/D — 10 _1 
and t/D = 10" 1 . 



D/R t < and anti-clockwise if D/R t > 0. Angles for 
which swimmers have linear trajectories, corresponding 
to a vanishing angular displacement in Fig. E{a) , are indi- 
cated by arrows. The purely rotating swimmer, with van- 
ishing \Ar\ but finite AO is indicated by the asymptote 
located at p « 103°. 

Apart from giving a useful insight into the circle swim- 
ming exhibited by the model, Fig. [5] can be used as a 
read-out to choose a particular swimmer, depending on 
the desired radius of trajectory and sense of motion. This 
can be useful given the wide range of radii that the tra- 
jectories can adopt. While the range in /? that results in 
clockwise motion is rather narrow, one can always obtain 
the desired direction of motion by considering negative 
values of /3 (mirror image swimmer). 



4 Velocity fields 

Based on the qualitative change in the trajectory with the 
swimmer angle, we expect that the flow field created by 
the swimmers also changes with /?. This is interesting in 
terms of the hydrodynamic interactions with surrounding 
objects, e.g. passive particles, both small tracers and ex- 
tended objects, and other swimmers. The average velocity 
field is calculated as 



f 

T 



Fjdt, 



(15) 



the swimmer with the links oriented at an angle /3/2 with 
respect to the x-axis. 

Figs. [H] and [7] depict the direction (arrows) and mag- 
nitude (colours) of v, for /3 = {0°, 53°, 107°} and /3 = 
{113°, 120°, 140°}, respectively. We plot the velocity field 
at two different length scales: on the left we show the range 
— 2D < x,y < 2D, and on the right we show the range 
— 20D < x, y < 20D. At small length scales the strength 
of the velocity field is higher close to beads 1 and 3, and is 
weakly dependent on j3. The flow pattern resembles that of 
the linear swimmer for small /3, for f3 > 107° the swimmers 
exhibit a recirculation pattern, up to distances comparable 
to the swimmers body length. This behaviour, however, is 
rapidly lost at larger length scales, as shown in the left 
panels of Figs. |5] and [7] 

We consider now the far-field behaviour. Quadrupo- 
lar decays (a velocity field which scales with distance as 
r~ 3 ) arise generally for swimmers whose stroke is invariant 
under a combined time-reversal and parity transforma- 
tion 38 . This invariance must be reflected in the velocity 
field - resulting in odd-power decays. For our swimmer the 
symmetry holds when the extension and deformation of 
both legs are identical. While the flow patterns for /3 > 
are in general reminiscent of a stresslet velocity field, a 
closer inspection of the magnitude of the velocity shows 
that it does indeed have a power-law decrease with dis- 
tance governed by an exponent n = —3. Fig. [SJa) plots 
the apparent decay exponent, n, measured at a large dis- 
tance from the swimmer, as a function of /?. We consider 
the decay along the the x- and y-axes, and define n as 



n(y = 0) 



d In I 



dh\x 



|j/=0) 



and 



d In I v . I 
n(x = 0) = — — % - x=0 
amy 



(16) 



(17) 



where i denotes a point in space with position vector 
in the frame of reference of the hydrodynamic centre of 



Since the exponent is measured at a finite distance from 
the swimmer, its value is in general non-integer. However, 
for sufficiently long measuring distances we expect to re- 
cover a single integer value for n. Along the y-axis the 
apparent exponent is always closer to n — —3, indicating 
that the quadrupolar term in the velocity field dominates 
at long distances as expected. This holds along the x-axis 
as well, except for f3 w {83°, 126°}, where higher order 
terms in the expansion dominate the decay, as indicated 
by the larger exponent n < —4. This behaviour is interest- 
ing, as it suggests that the symmetry of the swimmer at 
these angles suppresses the contribution of the quadrupo- 
lar term and gives way to higher order terms in the mul- 
tipole expansion. 

In order to demonstrate the change in the behaviour 
of the far field velocity as the swimmer loses its invariance 
under a time-reveral and partity transformation we have 
measured the apparent exponent for asymmetric swim- 
mers, where l\ and 1% have maximum extensions 0.8D 
and D, while all other parameters are kept as before. The 
results, presented in Fig. E£b), show apparent exponents 
along the x- and y-axes, which are consistent with a dipo- 
lar decay (n — —2). 
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5 Conclusions 

In this paper we have proposed a simple three-sphere model 
for a circle swimmer which is a natural generalisation of 
the linear Najafi-Golestanian swimmer [5]. The spheres 
are placed on a triangle such that two subsequent strokes 
are performed at an angle. We find that the radius and the 
sense of rotation of the swimmer trajectory depend del- 
icately on the separation angle between the rods joining 
the beads. The velocity field produced by this simple cir- 
cle swimmer exhibits a characteristic inverse-power decay 
at large distances. For swimming strokes invariant under 
a combined time-reversal and parity transformation we 
recover the expected quadrupolar decay for the velocity 
field, except for a narrow range of separation angles where 
a stronger decay is observed. As expected, for asymmetric 
swimming strokes, which do not posses the time-reversal 



and parity symmetry, we recover a decay consistent with 
a dipolar velocity field. 

The model can serve as a simple reference to help un- 
derstand more complex situations such as a circle swim- 
mer in a confining geometry or the collective properties of 
many circle swimmers. A further extension of the model 
is to swimmers whose angle varies in time [ZE]- Here we 
expect that more complicated modes of motion arise, de- 
pending on the interplay between the timescale of the link 
variation and the timescale of the swimmer angle varia- 
tion. Since the latter essentially controls the radius of the 
swimming trajectory, a slow variation of the swimmer an- 
gle relative to the links result in more complex swimming 
patterns, as shown in Fig. |Hl where the trajectories have a 
varying curvature leading to meandering motion. Any ar- 
bitrary two-dimensional trajectory can be generated by a 
suitable choice of P(t). Hence /3(i) can be used as a 'steer- 
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ing wheel' to navigate at will. Similar ideas, related to the 
controllability of the swimmer trajectories, have been re- 
cently explored by Alouges et al. [H] . The motion of circle 
swimmers under shear (e.g., linear shear or Poiseuille flow) 
can lead to new trajectories such as cycloids. Other exter- 
nal driving that could alter trajectories include a gravita- 
tional field or magnetic or electric fields [43]. One is just 
beginning to understand the topologies of these trajecto- 
ries for simple noise- free circle swimmers |36| . 

The trajectories of the circle swimmers studied in this 
paper could also be verified experimentally, e.g., by con- 
trolling the trajectories of immersed microbeads using op- 
tical traps, as performed by Leoni et al. |45| for the linear 
three-bead swimmer, and there is no obstacle in principle 
to do this for our circle swimmer as well. 

Finally, putting together more than three spherical 
beads could lead to more complicated motion, in three di- 
mensions, than the simple circular trajectories discussed 
in this paper. This results from the intricate translation- 
rotation coupling for biaxial particles. It was shown that 
the simple Brownian circle swimmer in three spatial di- 
mensions possesses a wealth of different trajectories with 
the circular helix being the simplest one [46 . It would be 
interesting to generalise our model to four spheres which 
are not in a common plane in order to access these com- 
plicated types of motion. 
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